Can spatial sorting associated with spawning migration explain evolution of body size and vertebral number in Anguilla eels?

Spatial sorting is a process that can contribute to microevolutionary change by assembling phenotypes through space, owing to nonrandom dispersal. Here we first build upon and develop the “neutral” version of the spatial sorting hypothesis by arguing that in systems that are not characterized by repeated range expansions, the evolutionary effects of variation in dispersal capacity and assortative mating might not be independent of but interact with natural selection. In addition to generating assortative mating, variation in dispersal capacity together with spatial and temporal variation in quality of spawning area is likely to influence both reproductive success and survival of spawning migrating individuals, and this will contribute to the evolution of dispersal‐enhancing traits. Next, we use a comparative approach to examine whether differences in spawning migration distance among 18 species of freshwater Anguilla eels have evolved in tandem with two dispersal‐favoring traits. In our analyses, we use information on spawning migration distance, body length, and vertebral number that was obtained from the literature, and a published whole mitochondrial DNA‐based phylogeny. Results from comparative analysis of independent contrasts showed that macroevolutionary shifts in body length throughout the phylogeny have been associated with concomitant shifts in spawning migration. Shifts in migration distance were not associated with shifts in number of vertebrae. These findings are consistent with the hypothesis that spatial sorting has contributed to the evolution of more elongated bodies in species with longer spawning migration distances, or resulted in evolution of longer migration distances in species with larger body size. This novel demonstration is important in that it expands the list of ecological settings and hierarchical levels of biological organization for which the spatial sorting hypothesis seems to have predictive power.     

RASPA: molecular simulation software for adsorption and diffusion in flexible nanoporous materials

A new software package, RASPA, for simulating adsorption and diffusion of molecules in flexible nanoporous materials is presented. The code implements the latest state-of-the-art algorithms for molecular dynamics and Monte Carlo (MC) in various ensembles including symplectic/measure-preserving integrators, Ewald summation, configurational-bias MC, continuous fractional component MC, reactive MC and Baker's minimisation. We show example applications of RASPA in computing coexistence properties, adsorption isotherms for single and multiple components, self- and collective diffusivities, reaction systems and visualisation. The software is released under the GNU General Public License.     

Chemotactic activity of pertussis toxin on murine lymphocytes. Its potential role on lymphocyte accumulation in the lung

Abstract The effects of pertussis toxin on lymphocyte migration were studied in vitro. In this study pertussis toxin significantly stimulated lymphocyte migration at concentrations of 0.1 and 1 μg ml⁻¹ using a microchamber and the leading-front method. Checkerboard analysis demonstrated that pertussis toxin causes directed migration of lymphocytes (chemotaxis). Heat-treatment of pertussis toxin abolished its capacity to cause this migration. When murine lymphocytes were preincubated with different concentrations of pertussis toxin, an inhibition of chemotaxis at the dosages of 0.1 and 1 μg ml⁻¹ was observed. On the other hand, lymphocytes derived from mice treated with pertussis toxin were not inhibited after subsequent exposure to pertussis toxin in vitro. Since lymphocyte accumulation in the lungs of mice treated with pertussis toxin has been well domenstrated, the results of our study could suggest a chemotactic activity of pertussis toxin in determining accumulation of lymphocytes in this organ.     

Migration timing and distance from shore of southbound eastern Pacific gray whales (Eschrichtius robustus) off Ensenada,Baja California,Mexico

Eastern Pacific gray whales were monitored off Ensenada, Mexico, during the southbound migration. The objectives were to determine southbound migration timing and width of the migration corridor during three seasons (2003–2006). Migration timing was determined by fitting a generalized additive model to the shore counts for each season and estimating the 10, 50, and 90 percentiles of the fitted curves. To estimate abundance from shore‐based counts, a probability density function for the shore based distances was estimated by a product of a gamma distribution fit to the boat survey distance data for 2006/2007 and a half‐normal detection function using combined data of the three seasons. The parameters of the gamma distribution were corrected to account for less boat survey effort carried out 20–40 km than 0–20 km from shore. The onset of the migration off Ensenada was in late December/early January and ended around 13 February. The median date was 23–26 January for the first and third season and a week early for the second season. Boat surveys indicated a wide (20 km) migration corridor but most gray whales traveled within 9.9 km from shore. The estimated total number of whales during watch hours was 2,298 (95% CI = 1,536–4,447).     

Force and torque modelling of drilling simulation for orthopaedic surgery

The advent of haptic simulation systems for orthopaedic surgery procedures has provided surgeons with an excellent tool for training and preoperative planning purposes. This is especially true for procedures involving the drilling of bone, which require a great amount of adroitness and experience due to difficulties arising from vibration and drill bit breakage. One of the potential difficulties with the drilling of bone is the lack of consistent material evacuation from the drill's flutes as the material tends to clog. This clogging leads to significant increases in force and torque experienced by the surgeon. Clogging was observed for feed rates greater than 0.5 mm/s and spindle speeds less than 2500 rpm. The drilling simulation systems that have been created to date do not address the issue of drill flute clogging. This paper presents force and torque prediction models that account for this phenomenon. The two coefficients of friction required by these models were determined via a set of calibration experiments. The accuracy of both models was evaluated by an additional set of validation experiments resulting in average R² regression correlation values of 0.9546 and 0.9209 for the force and torque prediction models, respectively. The resulting models can be adopted by haptic simulation systems to provide a more realistic tactile output.     

Exome resequencing reveals signatures of demographic and adaptive processes across the genome and range of black cottonwood (Populus trichocarpa)

Extant variation in temperate and boreal plant species has been influenced by both demographic histories associated with Pleistocene glacial cycles and adaptation to local climate. We used sequence capture to investigate the role of these neutral and adaptive processes in shaping diversity in black cottonwood (Populus trichocarpa). Nucleotide diversity and Tajima's D were lowest at replacement sites and highest at intergenic sites, while LD showed the opposite pattern. With samples grouped into three populations arrayed latitudinally, effective population size was highest in the north, followed by south and centre, and LD was highest in the south followed by the north and centre, suggesting a possible northern glacial refuge. F_ST outlier analysis revealed that promoter, 5′‐UTR and intronic sites were enriched for outliers compared with coding regions, while no outliers were found among intergenic sites. Codon usage bias was evident, and genes with synonymous outliers had 30% higher average expression compared with genes containing replacement outliers. These results suggest divergent selection related to regulation of gene expression is important to local adaptation in P. trichocarpa. Finally, within‐population selective sweeps were much more pronounced in the central population than in putative northern and southern refugia, which may reflect the different demographic histories of the populations and concomitant effects on signatures of genetic hitchhiking from standing variation.     

Novel Interaction of the Z-DNA Binding Domain of Human ADAR1 with the Oncogenic c-Myc Promoter G-Quadruplex

Both G-quadruplex and Z-DNA can be formed in G-rich and repetitive sequences on genome, and their formation and biological functions are controlled by specific proteins. Z-DNA binding proteins, such as human ADAR1, have a highly conserved Z-DNA binding domain having selective affinity to Z-DNA. Here, our study identifies the Z-DNA binding domain of human ADAR1 (hZα_ADAR1) as a novel G-quadruplex binding protein that recognizes c-myc promoter G-quadruplex formed in NHEIII₁ region and represses the gene expression. An electrophoretic migration shift assay shows the binding of hZα_ADAR1 to the intramolecular c-myc promoter G-quadruplex-forming DNA oligomer. To corroborate the binding of hZα_ADAR1 to the G-quadruplex, we conducted CD and NMR chemical shift perturbation analyses. CD results indicate that hZα_ADAR1 stabilizes the parallel-stranded conformation of the c-myc G-quadruplex. The NMR chemical shift perturbation data reveal that the G-quadruplex binding region in hZα_ADAR1 was almost identical with the Z-DNA binding region. Finally, promoter assay and Western blot analysis show that hZα_ADAR1 suppresses the c-myc expression promoted by NHEIII₁ region containing the G-quadruplex-forming sequence. This finding suggests a novel function of Z-DNA binding protein as a regulator of G-quadruplex-mediated gene expression.     

Wind and fairness in ski jumping: A computer modelling analysis

Wind is closely associated with the discussion of fairness in ski jumping. To counter-act its influence on the jump length, the International Ski Federation (FIS) has introduced a wind compensation approach. We applied three differently accurate computer models of the flight phase with wind (M1, M2, and M3) to study the jump length effects of various wind scenarios. The previously used model M1 is accurate for wind blowing in direction of the flight path, but inaccuracies are to be expected for wind directions deviating from the tangent to the flight path. M2 considers the change of airflow direction, but it does not consider the associated change in the angle of attack of the skis which additionally modifies drag and lift area time functions. M3 predicts the length effect for all wind directions within the plane of the flight trajectory without any mathematical simplification. Prediction errors of M3 are determined only by the quality of the input data: wind velocity, drag and lift area functions, take-off velocity, and weight. For comparing the three models, drag and lift area functions of an optimized reference jump were used. Results obtained with M2, which is much easier to handle than M3, did not deviate noticeably when compared to predictions of the reference model M3. Therefore, we suggest to use M2 in future applications. A comparison of M2 predictions with the FIS wind compensation system showed substantial discrepancies, for instance: in the first flight phase, tailwind can increase jump length, and headwind can decrease it; this is opposite of what had been anticipated before and is not considered in the current wind compensation system in ski jumping.